Fluctuations of the Fermi condensate in ideal gases 
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1. Introduction 

Starting quite a long time ago, Auluck and Kothari 1 , then May and finally 
Viefers, Ravndal, and Haugset j3j, discovered idependently that the specific heat of 
nonrelativistic ideal gases in two-dimensional (2D) boxes does not depend on the 
exclusion statistics. This interesting result eventually did not receive the attention it 
deserved until 1995, when Lee 0] rederived it by introducing an unified way of writing 
the thermodynamic properties of ideal gases in terms of polylogarithmic functions 5 . 
This formulation represented also an important extension of the Auluck and Kothari 
result and triggered further investigations (see for example Refs. jSlEIEI)- 

Since under canonical conditions a Bose and its corresponding Fermi gas 
are similar at the thermodynamic level, they have been called thermodynamically 
equivalent. If we denote by Cv (T, V, N) the heat capacity of a system at temperature 
T, volume V and particle number N, then the heat capacities Cv,i and Cy.2 of two 
thermodnamically equivalent systems are identical functions of T, V and N. Using 
this property, all the thermodynamic systems may be divided into equivalence classes 
|S] and by doing this one may observe that all the systems of ideal particles of the 
same constant density of states (DOS), but obeying Bose, Fermi, or even fractional 
exclusion statistics ^Oh belong to the same equivalence class [Jj]. 

The equivalence between Bose and Fermi gases was critically examined by Pathria 
[7|. He showed that the Lee's unified formulation of 2D ideal gases does not 
hold anymore below the Bose-Einstein condensation temperature of the Bose gas. 
Aparently, the 2D (or, more exactly, constant DOS) thermodynamic equivalence holds 
only above the Bose-Einstein condensation temperature. 

On the other hand, Crescimanno and Landsberg and one of us jH] showed 
that there is a one-to-one mapping between microscopic configurations of bosons, 
fcrmions or haldons, in systems with the same, constant DOS, which preserves the 
total excitation energy (i.e. the energy of the particles in the given configuration 
minus the energy of the system at zero kelvins is the same) . Based on this theorem, 
the thermodynamic equivalence of systems with equally spaced spectra should hold at 
any temperature in any detail, so Pathria's conclusion must be wrong. But what was 
overlooked there? 

The method of mapping microscopic configurations between systems of different 
exclusion statistics introduced in Ref. |5] for systems with constant DOS was 
extended in Ref. |12) for systems with any DOS and was called exclusion statistics 
transformation (EST). Systems connected by EST are thermodynamically equivalent 
by construction. If we take for example a Fermi system, transform it by EST into 
a Bose system, and then calculate independently the thermodynamical properties of 
these two systems by maximizing the entropy in each of them, for constant U (total 
energy) and N - i.e. assuming grandcanonical distribution on the single particle 
energy levels - we loose again the thermodynamic equivalence that we started with, 
at least is most of the cases ! The obvious conclusion is that one or both of these 
grandcanonical distributions lead to results in dezacord with the canonical ensemble. 
The question which of these grandcanonical distribution is closest to the canonical 
distribution is very dificult to answer, since ab initio canonical calculations are not 
easy to perform on general, large systems. 
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1.1. The Fermi condensate 

The concept of Fermi condensate was introduced in Refs. |121 ITS*] . For a general 
system, of spectrum consisting of single particle energy levels (we enumerate them 
such that eo < ei < . . .), one can calculate the grandcanonical probability, u>at , of 
having the lowest Nq energy levels occupied, the level No + 1 free and all the other 
energy levels with any occupation number |13| : 



Here Z is the partition function of the system and Z ex (N , (3, /3fi) = 
Tii^No+i {1 + ex P[ — P( e i ~ M)]} is the partition function of the levels iVo + 1, Ao + 2, . . .. 
Obviously, 



The probability distribution Q may have a maximum at, say A^max- The statistical 
interpretation of such a maximum is that in a physical system in contact with a 
particle reservoir, the lowest N (fv No,max) energy levels are always occupied. These 
No particles form the Fermi condensate. At any finite temperature, No is subject to 
fluctuations, denoted by 8 No- 

The configurations of fermions may be transformed by EST into configurations of 
bosons ^lEl- By this transformation the No degenerate fermions will be transformed 
into the Bose-Einstein condensate of the corresponding Bose system, and hence the 
name of Fermi condensate. For this reason the degenerate fermions will also be called 
Fermi condensate. For canonical Bose systems of constant DOS, the probability 
distribution of having No particles in the condensate, w N , have been studied in 
detail before (see for example Refill for a review). By construction, the canonical 
probability distribution of having Nq fermions in the condensate is also w^ Q ■ 

In Ref. ^31 it was shown that in a system of constant density of states, A^max 
- which corresponds to the maximum of wn _ is close to the average canonical 
ocupation of the ground state, but the two numbers are not the same. The distribution 
wn is not symmetric around the maximum, so (No) ^ A^max (where by (•) and 
(•) c we shall denote grandcanonical and canonical averages, respectively). Since if 
the Bose and Fermi condensates are separated from the rest of the particles, the 
grandcanonical Bose and Fermi distributions map onto each-other (see Section 5, Ref. 
|12|). the average number of particles in the Fermi condensate should be equal to the 
average ocupation of the bosonic ground state and in the Fermi system below the 
condensation temperature the usual Fermi distribution applies only to the particles 
above the condensate. For simplicity, to the fermions above the condensate we shall 
refer to as the particles in the thermally active layer. 

It is well known that Fermi canonical and grandcanonical ensembles are 
equivalent. In this paper we shall investigate the equivalence between the 
grandcanonical description of the Fermi gas and its correspondent canonical 
description in terms of the Bose gas. This is a new type of equivalence, first mentioned 
in Ref. ^2] an d which seems not to hold for general systems. Here we shall discuss 
the simplest systems, namely the ones with constant density of states, a, and we 
shall compare Wn and w N . An important parameter of the system is the quantity 
ak^T . For large ak^T we can do some analytical calculations, assuming that wn 




(1) 




(2) 
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Figure 1. Probability distribution of wpf as a function of jVrj. The maximum of 
the probability distribution is located at No trnax , which is given by the equation 
(fi — No,max/&)/kBT = Xo^aksT). For this particular plot ct/cbT = 100. 



has a gaussian shape, but, as we shall see in seection|21 this approximation is not too 
good for the evaluation of (S 2 Nq), the mean square fluctuation of Nq. To correct 
this and to extend our calculations to lower values of ak^T , in Section |3 we do 
numerical calculations of wm 01 w c N , (S 2 Nq) and (S 2 Nq) c . (For simplicity we use the 
somewhat simpler notations (S 2 N ) and (S 2 N ) C instead of the more rigorous ((5N ) 2 } 
and ((SN ) 2 ) C , respectively.) We reobtain the known result: (S 2 N ) C — ► ((2)(ak B T) 2 , 
i.e. the mean square fluctuation of No is of the same order as the number of particles 
in the thermally active layer or the number of particles on the excited states in the 
Bose gas. We also obtain that (5 2 Nq) — (S 2 Nq) c converge to a constant value, 0.39, as 
ak B T increases. This asymptotic behavior is proved analytically by the end of Section 

El 



2. Analytical evaluation of the fluctuations 



First we analyse Eq. analytically in the limit ak B T 3> 1. To do this we take 
logwjv , and transform all the summations into integrals. In this way we arrive to the 
expression [T3] 



logZ^ 



-0 



eo + A" (creo - 1) 



+ cr 



de log 



-/3(e- 



(3) 



where eo is the energy of the Nq single-particle level, the integral represents 
log Z qx (Nq, P, flu). The shape of the probability distribution is depicted in figure 
^ Since creo = N and d log Zn /cW = a^ 1 (d log ZN /deo), the value of e 
corresponding to the maximum of probability is given by the equation 

d\ogZ No 



de 



a I log 



f 



3 /9(co., 



c-/0 



(akBT)- 1 } 



= 0. 



(4) 
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and for ak B T > 1 Q3| 

eo,max = A* - k B Tlog[ak B T] . (5) 

We observe here that f3{/i — eo imax ) depends only on ak B T and not on /i, i.e. on the 
total particle number, as long as /i > k B T log[a ksT]. Therefore, as T decreases and /i 
becomes bigger than fc B T \og[ak B T] the probability distribution Q forms a maximum 
at Nq > 0. We say that at this temperature the condensate starts to form and the 
equation fi = k B T \og[ak B T] defines the condensation temperature. 

At low temperatures the maximum of wn becomes sharp and eo approaches /i. In 
this temperature range we shall approximate wn around the maximum by a gaussian 
distribution: 

{N - 7V ,max) 2 



w No w w(N 0jm3X ) ■ exp 



2A 2 



(6) 



The width of the gaussian is 
d 2 logZjs 



A 



'N 



dN 2 



No,, 



{akBT)- 1 - (1 + akBT)- 1 » {ak B T)- 2 . 



Equation © amounts to the use of the approximative function Z^ instead of Zn 



(7) 



z(a) 
-■No 



exp <{ log Z^ 



1 <9 2 logZ 



No X 2 - 



dN 2 



S 2 N } = Z No ^ x exp i - 



1 5 2 N 



To check the approximation we calculate first the total partition function @ as 

1 S 2 Nn 



Z^ = 



= Z 



d{SN )zP o a Z Nota 
nax v / 27rcrfc B T , 



2 (Cfc B r) 2 

m fun< 

d(5AT ) exp 



(8) 



2 (afc B T) 2 



(9) 



where 

log-Zjv 



= -f3[a 



-0,max 



^0,max 



l) + a 



/ de log 


l + e -^(e-M)' 


'SO, max 





akBT 



[(/3 M ) 2 - log 2 (afc B T)] - log^T) + afc B TL i2 {- exp[/3(/x - e ,max)]} , 

(10) 



where Li2 is Euler's dilogarithm function 5 . Using Eq. JSJ and the expansion 



i*a(2)L<_i 



(log 



6 (log|z|) 2 



(11) 



valid for any n > 1, we obtain the approximation 



log^ (a) 



crfc B r 



2 v ^ 6 
On the other hand, the exact partition function is 



/>oo 

/ de log 


1 + e -/3( 6 - M ) 
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which, if we apply again the approximation we get 



logZ = 



akBT 



y ( /3A1 ) 2 + _ (TfcB T. 



ak B T + log V2tt. (12) 
akBTLi 2 [- exp(j3fi)] , (13) 

(14) 
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The expansions of logZ and logZ^ are identical up to order ak B T '/ '(/3 fi) 2 , or 
ak B T / log 2 [Ck B T). The term logv^r ~ 0.92 from equation iff^l may be neglected, 
since is smaller than the order a k B T/ log 2 {ak B T) ^> 1, for ak B T 1. 

The fluctuation of the number of particles in the condensate, in the approximation 
®, is 

y/(8*N ) =ak B T. (15) 

On the other hand, the corresponding canonical fluctuation may be calculated for 
example by saddle point method applied to the equivalent Bose system [151 and 
gives 

V(5 2 N ) C = VW)ak B T , (16) 
where Q(x) is the Riemann Zeta function. Obviously, the two analytical 
approximations, 1|15|) and (|16|) . do not coincide and the question that remains to 
be answered is whether these distributions are indeed different, or simply the gaussian 
approximation (j^J is not good enough. 



3. Numerical evaluation of the fluctuations 



In this section we calculate the fluctuations numerically by introducing a recursion 
relation. From Eq. we obtain 

^jVo+i = exp[-/3(6Ar - /i)] 
w No 1 + exp[-(3{e No+1 - /j,)] 
and the value of iVo imax may be found by solving 
exp[-ff( £A r 0-max - fj,)] _ i 
1 + exp[-f3(e No max+ i - (j,)] 
If the density of states is constant and e^+i — £j = cr _1 for any i, then Eq. I|18|l becomes 
exp[-/3(e 0;max - fi)] = i ^ ig ^ 

1 + CXp[-/3(e ,max + 0- _1 - M)] 

Using Eqs. (|17l) and (|19fl we may now calculate numerically iV 0!max , (N ), and (5 2 N ). 
liak B T > 1, by writing exp[-/3(e 0!max +CT _1 -/i)] » exp[-/3(eo imax -^)](l-(fTfc B r) _1 ) 
Eq. (|19|) may be simplified to exp[/3(/x — eo,m a x)] = crk B T, which is the same as Eq. 
|J5J. Moreover, since around the maximum exp[/3(/i — eo, ma x)] 1, in the relevant 
energy interval we may transform Eq. i|17|) into 

^2±1 = {exp[/?(e A r - M )] + expHafceT)- 1 ]}- 1 « 1 - e^-M) + (^T)" 1 . (20) 

Let us now analyse the equivalent Bose gas. If again the system has a constant 
density of states we denote q = e- 1 /(' jk B T ) _ Then the canonical partition function for 
a system of N cx particles is JS] 



7 b 



In a canonical system of N particles, the probability w 1 ^ to have exactly iV ex particles 
in the excited states (not on the ground state) is proportional to Z\ — _ 1 |16| . 
so we have 

jv cx Wax -, 

fe=l y 
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Since N CK = N — N , the relative probability which corresponds to Eq. i(T7|) for 
fermions is 

w° 1 — a N ** 

= — = {l- exp[N m /(<Tk B T)]} exp^afceT)- 1 ] . (23) 

W N CX 1 

The most probable N cx is given by 

[1 - eM~PN c Ja)} exp[(afc B T)- 1 ] = 1 . (24) 

We want now to compare Eqs. (|17|l and ((23|) in the limit ak-^F ^S> 1. For this we 
take a Fermi and a Bose system with the same number of particles, AT. In the Fermi 
system we define ep = N/a (Fermi energy). We shall assume that both systems are 
below the condensation temperature and the number of particles in the condensate is 
No- Above the condensate we have N cx particles. For a condensed gas e-p — /.i < u^ 1 , 
so we can express N ex in Eq. I(23|) as iV ex = cr(ep — e ) = a(fi — e ). By doing so, Eq. 
lt2lH) becomes 



' M N + 1 



{1 - exp[/3(e„ - fi)]} exp^T)- 1 ] 



m 1 - exp[ ( 5(e„ - fj)] + (akBT)- 1 , (25) 

which is identical to Eq. (|20|l . Therefore the two probability distributions wn„ and 
^iVo a PP roacn each-other in the limit of large systems, i.e. when ak^F ^> 1. 

The numerical calculations, based on Eqs. (|17|) and 1)23(1 are plotted in figure [21 
We can observe that already for uk^F bigger than 1, the fluctuation of the particle 
number in the condensate is almost the same for the canonical Bose and grandcanonical 
Fermi systems. This justify the approach taken in Ref. ^3]> and for uk^F ^> 1, No 
may be calculated directly as the average number of particles in the Bose condensate, 
rather than by Eq. Q. 

Noticeble relative differences between canonical and grandcanonical results appear 
only for ak^F about 1 or below. For these values of ak^F the fluctuations depend 
on the exact location of /i, with respect to the single particle levels. For example let's 
say that [i £ (ejv-i, ejv), where ejv-i and en are two consecutive energy levels. In the 
limit — ejv-i) — > oo and for No = N, equation 1)17(1 becomes 

« exp[/?0u - ejv -i)] . (26) 

WN-l 

For /?(// — ejv-i) —> oo we can calculate (No) and (5 2 Nq) by taking into account only 
the levels e^-i and e^v and we obtain 

_ expjgg - «„_,)!(* 4-1) + AT _ N+1 _„ pl _^ tN _ im) 

exp[P(p - £Ar_i)J + 1 

and 

2 ^ = exp[/3(/i - ejy-i)] exp[-2/3(/x - ejy-i)] + (1 - exp[-/3(/x - fjy^i)]) 2 
^ ° ; exp^C/i - e N -i)] + 1 

1 - exp[-/3(/x - ejy-i)] 



exp[/3(/x - ejv-i)] + 1 
Therefore, for any \i 6 (e^r, e^v+i), 



exp[-/3(/i-eiv-i)] (28) 



hin jViVo) 1 / 2 [c 1/2 (2)afc B T] j=0- (29) 
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Figure 2. Numerical calculations of (<5 2 JV} 1/2 /[C 1/2 (2)°-A;bT] (a) and (S 2 N) 1 / 2 
(b). For low values of o-fcgT, (6 2 N) depends on the exact location of fi, with 
respect to the single particle levels. If ejv and ejv+i are two consecutive energy 
levels (ejv+i — ejv = CT_ 1 )i then the solid curves in both (a) and (b) figures 
correspond to fi = ejv + 0.1 ■ icr~ 1 (i = 0, 1, . . . , 5) from top to bottom. The 
dashed lines represent the canonical fluctuations, and the straight lines (y = 1 in 
a and y = £ 1//2 (2)x in b) correspond to the asymptotic values. 



The situation is different if for example fi = ejv- Then, applying the same 
algorithm as above, we get (N ) = N + 0.5 and (d 2 N)^ 2 = 0.5. In fi gure 2 we plotted 
(S 2 N)^ 2 /[C 1/2 {2)ak B T} (a) and {5 2 N) 1 / 2 (b) for M = e N + 0.1 ■ ia' 1 (t = 0, 1, . . . , 5). 
The fluctuations normalized to the asymptotic value, C}' 2 (2)akBT, are quite different 
for cfceT < 1, but the absolute values of the fluctuations are very close for any ak^T 
for both types of systems and any choice of p.. 

We notice also in figure [21(b) that although the difference -J (5 2 No) — J (5 2 Nq) c 
is very small for any ak^T ', it does not converge to zero as ak^T — > oo. Numerically 
we obtain 



y/(8 2 N a ) - ^(5 2 N a ) c w 0.39 for ak B t > 1 . (30) 

To explain this difference, let us note that the fluctuation of Nq in the grandcanonical 
ensemble, <5-/Vo = A^o — (No), may be viewed as the superposition of fluctuations 
comming from two sources: the canonical fluctuation of Nq around its average value, 
corresponding to the total particle number N, denoted by S N Nq(= No — (No)n), and 
the fluctuation of (No) n due to the grandcanonical fluctuation of N. Assuming small 
fluctuations, the variation of (No) jv due to the variation of N may be written as 

s ( n o)n = dN -SN. 
Collecting all these together we write 

SN =N - (N ) = N - (N a ) N + (N )n - (N ) - (W C + ^J^^ N ■ (31) 

Below the condensation, d(No)/dN — 1 (temperature stays constant). Moreover, 
well below the condensation temperature SNq and 5N are independent, since the 
condensate may be viewed as a reservoir of particles of zero energy |17j , and Eq. (|31[) 
leads to 

(6 2 N ) = (8 2 N ) C + (S 2 N) . (32) 
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For high enough ak B T and fin we have (S 2 N) « crk B T, which, if plugged in Eq. 
gives 

(S 2 N ) « VC(2)(afcBT)2 + d B T « VW)^ B T + — -L= . (33) 

As we expect, (2 V /C(2)) _1 ~ 0.39. 

Using the same method we calculate the fluctuation of the number of particles 
in the thermally active layer without doing any extra numerics. Again we denote by 
(N ex ) the average number of particles in the thermally active layer and the fluctuation 
SN ex can again be written as 

SN CX = N cx - (N cx ) = N cx - (N CX ) N + (N CX ) N - (N cx ) = 5 C N CX + ^^-SN . (34) 

By (N ex ) at we denote the average number of particles in the thermally active layer at 
fixed N. Well below the condensation temperature (N ex )jsr does not depend on TV, so 
from Eq. |53J we get SN CX = S C N CX and 

(5 2 N cx ) = (5 2 N cx ) c = (6 2 N a ) c . (35) 



4. Conclusions 



The thermodynamic equivalence between ideal bosons and fermions with the same 
constant density of states a ^ |21 El ED E3 EH E3 apparently is lost below the Bose- 
Einstein condensation temperature [?]• On the other hand it was proven that if the 
Bose and the Fermi systems have the same spectrum consisting of nondegenerate, 
equidistant single particle states (like for example particles in a one-dimensional 
harmonic potential), then the canonical thermodynamic equivalence between the two 
systems is preserved down to zero temperature in the smallest details m Ej • This 
apparent contradiction is due to the fact that below the condensation temperature T c , 
in the Fermi gas the values of intensive parameters like the chemical potential and also 
the population of single particle levels in the canonical ensemble arc changed slightly 
from their corresponding values in the grandcanonical ensemble. 

Below T c , to the Bose-Einstein condensate in the Bose system it corresponds in 
the Fermi system a degenerate gas of the same number of particles, at the lowest part 
of the spectrum. Because of this correspondence the degenerate fermionic subsystem 
is called here the Fermi condensate The Fermi condensate is manifested also 

in the probability distribution of the grandcanonical ensemble. One can calculate 
the probability to have N degenerate particles (see Eq. 1). Below T c this probability 
distribution has a maximum for Ao > and we showed numerically and analytically in 
Section |3 that the grandcanonical average of No is the same as the canonical average. 

In this paper we did both, analytical and numerical calculations of the number 
of particles in the condensate and in the thermal active layer. We calculated also 
their fluctuations. The grandcanonical fluctuation of No is almost the same as the 
fluctuation in the canonical ensemble. Although the average values (No) and (iVo) 
are identical, for large values of ak B T the fluctuations (5 2 No) 1 ^ 2 and ((S 2 No) c ) 1 ^ 2 
differ by a small, but constant value, 0.39. This is due to the extra contribution to 
(8 2 No) 1 ^ 2 1 given by the grandcanonical fluctuation of the total particle number. 

The fermions in the thermally active layer correspond to the bosons on the 
excited energy levels. Canonical and grandcanonical averages of N cx are the 
same. Moreover, well below the condensation temperature, where Maxwell Deamon's 
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ensemble is applicable J7|, the fluctuation of iV ex is the same in both, canonical and 
grandcanonical ensembles i|35[l . 
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